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Charmed meson decays to a light pseudoscalar (P) and light vector (V) meson 

are analyzed taking account of 77-77' mixing. A frequently-used octet-singlet 

mixing angle of 19.5° is compared with a value of 11.7° favored by a recent 

analysis of D — > PP decays. 
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Decays of the charmed mesons D , D + , and Df to a light pseudoscalar meson P and a 
light pseudoscalar meson V were analyzed within the framework of flavor SU(3) in Refs. [1] 
and [2]. A frequently-used octet-singlet mixing angle between 77 and 77' of 9 V = 19.5° was 
CSj ! used in Ref. [1], while Ref. [2] used 9 V = 14.4° based on a recent KLOE analysis [3]. In a 

study of -D( s ) — > PP [4], a best fit to Cab ibbo- favored decay rates was found for 9 V = 11.7°. 
^ I In the present Brief Report we update fits to D( s ) — ¥ PV decays including two decay modes 

not considered in [1], and compare fits based on 9 V = 19.5° and 11.7°. 

We review our notation [4]. The angle 9 V describing octet-singlet mixing between 77 and 
77' is defined by 



(N 



in 

o 
o 



77 = — 778 cos Or, — 771 sin 9 V , T]' = —r)gsin9 ri + r)iCOs9 TI , where (1) 

77s = (uu + aid — 2ss)/v6 , 771 = (uu + dd + ss)/v3 , (2) 

Our previous analysis of PV decays utilized 9 V = arcsin(l/3) = 19.5°, for which 

77 = (ss — uu — dd)/y3 , n' = (2ss + uu + dd)/V6 . (3) 

We consider also 9 V = 11.7°, for which an exact fit was found in Ref. [4] to Cab ibbo- favored 
decays. 

We refer to Ref. [1] for notation. Amplitudes defined there include color-favored tree 
(T), color-suppressed tree (C), exchange (E), and annihilation (A), with a subscript P or 
V denoting the meson containing the spectator quark. Fitting the Cabibbo-favored data 
quoted there, we found two solutions ("A" and "B"), distinguished by \Ty\ < \Cp\ (A) and 
\Ty\ > \Cp\ (B). Fits to singly- Cabibbo-suppressed data then favored solutions consistent 
with the set "A," which we shall consider from now on. In Table I we show the results of 
this fit. 

We then fit amplitudes involving 77 and 77', obtaining values for the amplitudes Tp, Cy, 
and Ey. These are compared for the two- most-favored solutions (denoted by Al and A2) 
in Tables II and III. Predictions for the branching fraction B(D° — > K*°i]'), listed in the 
last columns of Tables II and III, differ slightly between solutions Al and A2. 



Table I: Solution in Cabibbo-favored charmed meson decays to PV final states favored by 
fits [1] to singly- Cabibbo-favored decays. 



PV 
amplitude 



Magnitude 
flO" 6 ) 



Relative 

strong phase 



C P 
Ep 



3.95 ±0.07 
4.88 ±0.15 
2.94 ±0.09 



S, 



CpTy 



5 



EpT v 



-162 ±1)° 
-93 ±3)° 



Table II: Solutions for Tp, Cy, and Ey amplitudes in Cabibbo-favored charmed meson 
decays to PV final states. Solutions Al and A2 correspond to \Ty\ < \Cp\. Here the 77 — 77' 
mixing angle is 9 V = 19.5°. 



No. 


PV 


Magnitude 


Relative 


B(D° ->■ K*°r]') 




ampl. 


(10" 6 ) 


strong phase 


(10- 4 ) 


Al 


T P 


7.46±0.21 


Assumed 






Cy 


3.46±0.18 


5 Cv t v = (172 ±3)° 


1.52 ±0.22 




Ey 


2.37±0.19 


<Wv = (-H0±4) 




A2 


T P 


6.51±0.23 


Assumed 






Cy 


2.47±0.22 


5evr P = (-174±4)° 


1.96 ±0.23 




Ey 


3.39±0.16 


Se v t v = (-96 ± 3)° 





Here we use a line of thought different from the analysis of Ref. [1]. We now calculate 
global x 2 values for fits to singly- Cabibbo-suppressed D° — > PV decays for the solutions 
Al and A2. We compare the x 2 values for the fit with 8 V = 19.5° [including branching 
fractions B(D° ->■ rju) = (0.221 ± 0.023)% [6] and B(D° ->■ t?0) = (1.4 ± 0.5) x 10~ 4 [5] 
omitted in the original article] with values for a fit to the same data with 6 V = 11.7°. These 
results are shown in Tables IV and V, respectively. 



We see that solution Al is favored for both 6 r , 



19.5° and 0„ 



11.7°. The solution 



A2 is disfavored since its prediction for B(D° — > i]<j)) is much higher than the experimental 
value in both cases. The same conclusion is reached in Ref. [2] for 6 V = 14.4°. We will now 
disregard the A2 solution and only use the Al solution for the rest of the analysis. 

The next step is to use observed Cabibbo-favored decays to obtain the annihilation 
amplitudes Ap and Ay using the amplitudes for D s — > (K* K + ,K K* + ,n + uj) (as quoted 
in Table VI) and the Al solutions. Since we use only 3 independent inputs to obtain 4 
independent parameters (real and imaginary parts of Ap and Ay) instead of obtaining 
unique solutions, we obtain a zone of allowed parameter space. We first form a grid of 
\A P \ and \A V \ values, and for every point on this grid, use the amplitudes for D s — > 
(K*°K + ,K°K* + ) to obtain the phases of Ap and Ay relative to Ty (assumed real, as 
previously.) Thus for every point on this grid we now have an amplitude for the decay 
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Table III: Same as Table II but with 6 V = 11.7°. 



No. 


PV 


Magnitude 


Relative 


B(D U -»■ K*°r]') 




ampl. 


(1(T 6 ) 


strong phase 


(io- 4 ) 


Al 


T P 


7.69±0.21 


Assumed 






C v 


4.05±0.17 


5 CvTv = (162±4)° 


1.19 ±0.12 




Ey 


1.11±0.22 


5 EvTv = (-130 ± 10)° 




A2 


T P 


5.68±0.23 


Assumed 






C v 


1.74±0.23 


5 Cv t p = (-162 ± 6)° 


2.19 ±0.16 




Ey 


3.82±0.15 


Sc v t v = (-87 ±3) c 





Table IV: Global x 2 values for fits to singly-Cabibbo-suppressed D° — > PV decays. Also 
included are the process that contribute the most to a high x 2 value. Here we have taken 



9 V = 19.5°. 



No. 


Global 




Worst Processes (Highest Ax 2 values) 






x 2 


Decay Channel 


B th (%) B expt (%) 


A X 2 


Al 


55.9 


D° -*7]<f) 




(4.0 ± 0.4) x 10~ 2 (1.4 ± 0.5) x 10" 2 


16.8 






D° — > T)U> 




0.33 ±0.02 0.221 ±0.023 


11.3 


A2 


82.4 


D° -+Tj<j) 




(5.9 ± 0.4) x 10~ 2 (1.4 ± 0.5) x 10~ 2 


45.8 






D ^n°p 




0.27 ±0.02 0.373 ±0.022 


10.1 



D s — > ir + u using the amplitude representation from Table VI. We now select only those 
points on this grid that are allowed by the experimental value \A(D S — > n + u)\ including 
its one-sigma error bar. 

Since there is a two-fold discrete ambiguity in choosing the phase of Ap relative to Cy 
or that of Ay relative to Cp, we obtain four different sets of allowed zones on the parameter 
space defined by \Ap\ and \Ay\. The allowed zones for T) = 19.5° and 9 V = 11.7° are shown 
in Fig. 1 and Fig. 2 respectively. One may associate unique phases with Ap and Ay (that 
may be determined following the method explained above) for every point on the \Ap\ 
\Ay\ plane in each of these figures. 



Table V: Same as Table IV but with 9 V = 11.7°. 



No. 


Global 


Worst Processes (Highest Ax 2 values) 






x 2 


Decay Channel Bth(%) B cxpt (%) 


A X 2 


Al 


35.8 


D°^n + p- 0.39 ±0.03 0.497 ±0.023 


8.5 






D° ^ricu 0.30 ± 0.02 0.221 ± 0.023 


6.7 


A2 


131.4 


D°^r](j) (9.2 ±0.6) x 10~ 2 (1.4 ±0.5) x 10~ 2 


100.0 






D°^7r°p 0.27 ±0.02 0.373 ± 0.022 


11.4 



Table VI: Branching ratios [5] and invariant amplitudes for Cabibbo-favored decays of D s 
used to obtain Ap and Ay. V is the 77 — 77' mixing angle. <\>\ = arcsin(l/v^3) = 35.3°. 



Meson Decay 


Representation 


B[5] 


p* 


1*41 


mode 




(%) 


(MeV) 


(10- 6 ) 


D+ K* U K + 


Cp + A v 


3.9 ±0.6 


682.4 


3.97 ±0.31 


K°K*+ 


Cy + Ap 


5.3 ±1.2 


683.2 


4.61 ±0.52 


7l + UJ 


j= 2 (A v + A P ) 


0.25 ±0.09 


821.8 


0.76 ±0.14 


P + V 


T P cos(fl T; + 00 - ^^ sin(0„ + 00 


13.0 ±2.2 


723.8 


6.63 ±0.56 


p + T)' 


T P sin(0„ + 00 + ^±Av cos(0„ + 0i) 


12.2 ±2.0 


464.8 


12.5 ±1.0 



Table VII: Range of predicted branching ratios for D s — > (77,77') p + using both $ v = 19. 5 C 



and 9 V - 


= 11.7°. 














Decay 
mode 


9 V = 
Min(%) 


19.5° 
Max(%) 


Or) = 

Min(%) 


11.7° 
Max(%) 




B(D S -+ rjp + ) 
B(D S -v v 'p+) 


3.80 
2.71 


6.39 

3.41 


6.27 

2.45 


8.35 
3.04 



1.16 c _ /_-. p-i+83^0 
1.04? °A P — { 1 '- )1 -74J 


(4) 


0.34 e / iq+10\o 
0.31? °A V - (-19_9 ) 


(5) 


■■ (5.6 ± 1.2)% and B(D S ->■ r]' p+) 


= (2.9 ± 1.2)%. 



To conclude one may now use the range of possible Ap and Ay values to predict B(D S — > 
r]p + ) and B{D S — > r]'p + ). In Ref. [1] we used the solution: 

\A P \ = 1.36J 
\A V \ = 1.25" 

which led us to obtain B(D S — > r]p + ) = 

Over the region of allowed values for Apry), the central values for these Cabibbo-favored 

D s branching ratios vary over the ranges shown in Table VII. 

The predictions for B(D S — > i]p + ) are a bit higher in the new fit using 9 V = 11.7° and 
slightly closer to the experimental value [5] quoted in Table VI. No improvement is seen 
in the prediction for B(D S — ¥ i]'p + ) in the new fit using 9 r) = 11.7°. The experimental 
values for both these branching ratios [5], as quoted in Table VI, are much higher than the 
predictions using this analysis. As mentioned in Ref. [1], the relation 

\A(D S -► pV)| 2 = \T P \ 2 + \A(D S -». n + co)\ 2 - \A{D S -»- p + r 1 )\ 2 (6) 

is very badly obeyed with the present values of B(D S — > i]p + ) and B(D S — >■ r]'p + ), leading us 
to suspect either that they have been overestimated experimentally, or that disconnected 
diagrams (as studied in [4]) play a larger role than anticipated. The scarcity of available 
data for Cabibbo-favored processes prevents such an analysis in the D — > PV case. 

This work was supported in part by the United States Department of Energy under 
Grant No. DE-FG02-90ER40560. 
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Figure 1: Allowed values for \Ap\ and \Ay\ for V = 19.5°. In order to obtain the phase of 
Ap(y) we either add (denoted by +) or subtract (denoted by — ) its phase relative to Cv(p) 
from the phase of C\/(p)- Clockwise from top left the 4 panels represent the phase choices: 

a) ++, b) H — , c) and d) — h Thus one may associate a unique value of the phase of 

Ap(y) with every point on the parameter space in these plots. 
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Figure 2: Same as Fig. 1 but with 9„ = 11.7°. 
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